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ABSTRACT 

We present a first measurement of the stellar mass component of galaxy clusters selected via 
the Sunyaev-Zel'dovich (SZ) effect, using 3.6 [J.m and 4.5 \xm photometry from the Spitzer 
Space Telescope. Our sample consists of 14 clusters detected by the Atacama Cosmology 
Telescope (ACT), which span the redshift range 0.27 < z < 1.07 (median z = 0.50), and 
have dynamical mass measurements, accurate to about 30 per cent, with median M500 = 
6.9 x 10 14 Mq. We measure the 3.6 [im and 4.5 |xm galaxy luminosity functions, finding the 
characteristic magnitude (m*) and faint-end slope (a) to be similar to those for IR-selected 
cluster samples. We perform the first measurements of the scaling of SZ-observables (F500 and 
yo) with both brightest cluster galaxy (BCG) stellar mass and total cluster stellar mass (Mf^). 
We find a significant correlation between BCG stellar mass and I500 (E{z)~ 2 ^ D^Y^qq °c 
Ml' 2±Q ' 6 ), although we are not able to obtain a strong constraint on the slope of the rela- 
tion due to the small sample size. Additionally, we obtain E(zY 2 ! 3 D 2 A Y 5m oc Mf™ 10±a6 for 
the scaling with total stellar mass. The mass fraction in stars spans the range 0.006-0.034, 
with the second ranked cluster in terms of dynamical mass (ACT-CL J0237-4939) having 
an unusually low total stellar mass and the lowest stellar mass fraction. For the five clusters 
with gas mass measurements available in the literature, we see no evidence for a shortfall of 
baryons relative to the cosmic mean value. 

Key words: cosmology: observations - galaxies: clusters: general - galaxies: luminosity 
function, mass function - galaxies: stellar content 
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1 INTRODUCTION 

Galaxy clusters are the most massive gravitationally bound objects 
in the Universe, and as such their changing abundance with redshift 
traces the process of structure formation. By studying the evolution 
of their properties with redshift, we can learn the assembly his- 
tory of their constituent dark matter, gas, and galaxies. The large 
masses of clusters ensure that the gas and stellar mass observed 
within them remain gravitationally bound. They therefore represent 
a fair sample of the Universe as it evolves over cosmic time. Mea- 
surements of cluster gas fractions have been used to constrain cos- 
mological parameters, including the dark energy equation of state 
(e.g., |Allen et al.|2008| >. While the hot gas in clusters makes up the 
majority of the baryonic component (around 80 per cent), accurate 
measurements of the mass fraction in the minority stellar compo- 
nent are also required to gain insight into the physical processes 
occurring within clusters. 

It is known that the physics which determines the properties of 
the intracluster medium (ICM) is more complicated than simply the 
action of gravitational collapse alone, as the observed scaling rela- 
tions, such as between X-ray luminosity and temperature, deviate 
from the self-similar expectation (e.g., Markevitch 1998 , Arnaud & 



Evrard 1999; Pratt et al. 2009), indicating that an additional source 
of energy is heating the ICM. While some energy is injected by 
supernovae (SNe) within galaxies, it is likely that the bulk of the 
energy comes from Active Galactic Nuclei (AGN) in the centres 
of clusters, as observations of low redshift clusters show that AGN 
jets, seen in radio imaging, carve out cavities in the hot gas ob- 
served at X-ray wavelengths (e.g., |Birzan et al.]2 004; Mc Namara| 
|et al.|2005|rBlanton et al.l20TT) . Feedback by AGN and SNe must 
also leave an imprint on the galaxy populations of clusters, and is a 
crucial component of galaxy formation models (e.g., |De Lucia et aT] 
|2006| [Bower et a l. 2006, 2008), where it is invoked to quench star 
formation in massive haloes. Therefore, in principle, it is possible 
to constrain the strength of feedback processes by measuring the 
stellar fractions of clusters (e.g., Bode, Ostriker & Vikhlinin 2009 ). 

The baryon fractions of clusters in the local Universe have 
been measured by several studies (e.g., Lin et al. 2003 Gonza-| 
|lez et al.|2007||Andreon|20T0l|Balogh et al.|2011| > They find that 
the mass fraction contained in stars is smaller for more massive 
clusters, while conversely the gas fractions are larger (e.g.. Gon- 
|zalez et aT| [2007). Within Rsqq (the radius at which the enclosed 
density is 500 times the critical density of the Universe), there 
may be a shortfall of baryons with respect to the cosmic value in- 
ferred from measurements of the cosmic microwave background 
(/b = 0. 166± 0.020;|Komatsu et al.|2011), with cluster studies find- 
ing /b ~ 0.13 (e.g., |Gonzalez et al.|2007| >, although the uncertain- 
ties are large, and subject to various systematic effects. There are 
few measurements at intermediate redshift: Giodini et al. ( 2009 1 de- 
scribe measurements of the baryon fractions in 0. 1 < z < 1 groups 
and clusters found in the COSMOS ( |Scoville et aT]|200 : 7) field, 
while [Lin et al.| ([2012l present estimates for a heterogeneous sam- 
ple of 45 clusters at 0.1 < z < 0.6 using archival X-ray data and 
infrared photometry from the Wide-Field Infrared Survey Explorer 
(WiSE; |WrTght et al.|2010) >. 

In this paper, we report a first characterisation of the stellar 
fractions, stellar mass scaling relations, and properties of brightest 
cluster galaxies (BCGs), using infrared data from the Spitzer Space 
Telescope, for a sample of galaxy clusters selected via the Sunyaev- 
Zel'dovich effect (SZ; Sun yaev & Zeldovich|1972| by the Atacama 
Cosmology Telescope (ACT; Swetz et al.|2011[ l The SZ effect is 
the inverse Compton scattering of cosmic microwave background 



photons by electrons in the hot gas atmospheres of clusters. In 
principle, it allows the construction of approximately mass-limited, 
redshift-independent samples of clusters (see, e.g., the review by 
ICarlstrom, Holder & Reese|2002| >. 

The structure of this paper is as follows. In Section [2] we 
briefly describe the cluster sample and associated ancillary data, 
which include dynamical mass estimates obtained from extensive 
spectroscopic observations with 8 m class telescopes (Sifo n et al.| 
|2012[ >. We also describe the processing of the Spitzer infrared imag- 
ing used in this work. In Section [3] we present a measurement of 
the 3.6 vim and 4.5 |i.m galaxy luminosity functions, and compare 
our findings with results obtained for IR-selected cluster samples. 
We examine the IR-properties of the BCGs, and their scaling rela- 
tions with dynamical mass and SZ observables, in Section [4] We 
similarly present scaling relations with total stellar mass, and the 
stellar fractions, in Section[5] We discuss our findings in Section|6] 
and conclude in Section|7] 

We assume a cosmology with Q. m = 0.3, Q.\ = 0.7, and 
Ho = 70 km s _1 Mpc~' throughout. All magnitudes are on the AB 



system ( |Oke|1974} , unless otherwise stated. We adopt a |Salpeterj 
{1955} Initial Mass Function (IMF) for stellar mass estimates. All 
estimates of cluster masses and SZ-signals are measured within a 
characteristic radius defined with respect to the critical density at 
the cluster redshift. 



2 DATA AND OBSERVATIONS 

2.1 Cluster sample and optical observations 



The cluster sample used in this work is drawn from the Marriage 
|et al.| j2011) SZ-selected cluster catalogue, detected in 148 GHz 
observations (conducted in 2008) with ACT. a 6 m telescope lo- 
cated in Northern Chile. The pipeline used to produce maps (di- 
rectly calibrated to WMAP; see |Hajian et al.|2011|> fro m the ACT 
timestream data, is described in |T5u^mer*et*aL 1 2013 ). Follow-up 
optical imaging observations of these clusters, using the NTT and 
SOAR telescopes, are reported in Menanteau et al. ( 2010|, and |Seh-| 
|gal et al.| l (201 ip discuss the cosmological constraints obtained from 
this sample. 

In this work, we use a subset of 14 clusters from the ACT sam- 
ple that have received extensive spectroscopic observations with 
Gemini South and the VLT ( |Sif6n et al.pOTI) . On average, red- 
shifts were measured for 60 galaxies per cluster, resulting in dy- 
namical mass estimates with a typical uncertainty of about 30 per 
cent. These measurements were used to investigate the mass scaling 
relation with the SZ-signal characterised using three different meth- 
ods. We use two of these types of measurements in this paper; both 
are described in detail in Hasselfield et al. (2013). The first is the 
dimensionless, matched filter SZ amplitude (yo), which is extracted 
from maps filtered using an Arnaud et al. (2010) profile with fixed 
angular scale 6500 = 5.9', and subsequently corrected to the angular 
scale corresponding to R500 at me observed cluster redshift as de- 
scribed in Hasselfield et al. Our second estimate is the widely used, 
spherically integrated Compton signal (I500), measured within the 
radius Rsoo- It is important to note that both of these measurements 
ultimately derive from the matched filter amplitude, and are there- 
fore highly correlated with each other and dependent upon the as- 
sumed |^maud^eFaL] |20 10) model for the cluster SZ signal. We 
use the updated I500 and yo values presented in the Appendix of 
|Hasselfield et al.| these were measured using improved ACT maps 
which include data obtained in the 2009-2010 observing seasons. 
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Figure 1. Redshift, dynamical mass, and 7500 distributions for the Sunyaev-Zel'dovich effect selected cluster sample use d in this work (see|Sifon et al.|2012| 
for details). The clear outlier in the Ysqq distribution is the z = 0.87 merger system ACT-CL J0102-4915 ("El Gordo"; see.Me nanteau et al.|2012a[ l 



Table 1. Properties of the cluster sample used in this work. Masses are dynamical estimates obtained from line of sight velocity dispersion measurements; all 
values have been rescaled from R200 ( as reported by Sifon et al. 2012 1 to R^oo- 



Name 


z 


M 50 o(10 14 M Q ) 


Ksooao- 11 ) 


yo (10- 4 ) 


S500 (Mpc) 


ACT-CL JO 102-49 15 


0.870 


9.8±2.3 


9.7±1.7 


7.17±0.88 


1.1 


ACT-CL J0215-5212 


0.480 


5.8± 1.7 


2.0±0.8 


1.18±0.27 


1.1 


ACT-CL J0232-5257 


0.556 


4.2±1.5 


1.0±0.6 


0.94 ±0.27 


0.9 


ACT-CL J0235-5121 


0.278 


7.0±2.1 


5.0±1.7 


1.14 ±0.22 


1.2 


ACT-CL J0237-4939 


0.334 


12.2±2.5 


2.6±1.8 


1.09 ±0.31 


1.4 


ACT-CL J0304-4921 


0.392 


7.7 ± 1.8 


5.8±2.2 


2.08 ±0.41 


1.2 


ACT-CL J0330-5227 


0.442 


10.3 ±2.4 


4.6±1.0 


1.68 ±0.24 


1.3 


ACT-CL J0346-5438 


0.530 


6.4 ±1.4 


2.8 ±0.9 


1.67 ±0.34 


1.1 


ACT-CL J0438-5419 


0.421 


12.7±3.0 


7.4±0.8 


2.07 ±0.17 


1.4 


ACT-CL J0509-5341 


0.461 


3.5± 1.4 


2.6±0.6 


1.13±0.19 


0.9 


ACT-CL J0528-5259 


0.768 


3.6±1.3 


0.8±0.3 


1.03 ±0.28 


0.8 


ACT-CL J0546-5345 


1.066 


4.8 ±2.5 


1.7±0.3 


2.54 ±0.39 


0.8 


ACT-CL J0559-5249 


0.609 


8.9±2.6 


2.4 ±0.5 


1.48 ±0.23 


1.2 


ACT-CL J06 16-5227 


0.684 


6.8 ±2.9 


2.5 ±0.5 


1.87 ±0.28 


1.0 



Note that because our IRAC data do not cover out to R200 
for the lower redshift objects in our sample (see Section [2~2| , we 
rescale the measurements of Sifon et al. (2012| > at i?200 to ^500 



Duffy et al.|p008) >. We 



using the concentration-mass relation of 
include the uncertainty introduced by the scatter in this relation in 
the error bars on the M500 measurements. This results in the frac- 
tional mass errors being approximately 9 per cent larger than the 
uncertainties quoted in Sifon et al. (2012). Table[T]lists the cluster 
properties. 



Fig. [T] shows the redshift, mass, and /500 distributions for the 
clusters used in this work. The sample spans the redshift range 
0.27 < z < 1.07, with median z = 0.50, and has median mass 
M500 = 6.9 x 10 14 Mq. The object with the largest intrinsic Y500, 
which is a clear outlier from the distribution shown in the right 
panel of Fig. [T] is the spectacular "El Gordo" (ACT-CL JO 102- 
4915) merger system at z = 0.87, for which we have previously 
published a detailed multiwavelength analysis iMenanteau et al.| 
|2012a| >. 



2.2 Spitzer IRAC imaging and photometry 

We obtained IRAC 3.6 p.m and 4.5 \im observations of the clus- 
ter sample during the period August-December 2010 (program ID: 
70149, PI:Menanteau), using a 2 x 2 grid of IRAC pointings centred 
on each cluster position. A total of 10 x 100 sec frames were ob- 
tained in each channel at each grid position, using a large scale cy- 
cling dither pattern. The basic calibrated data (BCD) images were 
corrected for pulldown using the software of Ashby & Hor^J and 
then mosaicked using MOPEX (Makovoz & Khan 2005) to give 
images which are ~ 13' on a side with 0.6" pixel scale. The mo- 
saics for each channel were then registered to a common pixel co- 
ordinate system. These mosaics cover out to R^qq for every cluster 
in the sample. By inserting synthetic point sources, we estimate the 
80% completeness depths of the final maps to be ~ 22.6 mag (AB) 
in both channels. 

Matched aperture photometry was performed on the IRAC 
maps using SEXTRACTOR (Berti n & Arnouts|1996| > in dual image 
mode, using the 3.6 (im channel as the detection band. We mea- 



1 See |http://irsa,ipac.caltech.edu/data/SPITZER/docs/| 
|dataanalysistools/tools/contributed/ irac/ f ixpulldown/| 
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Figure 2. Composite luminosity functions of ACT Sunyaev-Zel'dovich effect selected galaxy clusters in two redshift bins, for each IRAC channel. The solid 
curve shows the best fitting Schechter function with m* and a as free parameters, while the dashed curve shows the best fit with a fixed at —0.8. 



sure fluxes through 4" diameter circular apertures, which are cor- 
rected to estimates of total magnitude using aperture corrections (as 
measured by |Barmby et al.|2008) of -0.35 ± 0.04, -0.37 ± 0.04 
mag in the 3.6 [im, 4.5 \im channels respectively. The photomet- 
ric uncertainties were scaled upwards by factors of 2.8, 2.6 in 
the 3.6 (tm, 4.5 \im channels respectively, in order to account for 
noise correlation between pixels introduced in the production of 
the mosaics which is not taken into account in the SEXTRACTOR 
error estimates. These scaling factors were determined using the 
method outlined in Barmby et al. (2008). Finally, the uncertainties 
in the aperture corrections were added to the photometric errors in 
quadrature. 



3 INFRARED LUMINOSITY FUNCTIONS 

The luminosity function (LF) encodes key information about 
galaxy populations, and has been used to characterise the variation 
of galaxy properties with environment (e.g 



B lan ton et al. 



2003 



De Propris et al. 2003] |Croton et aL][2003| |Popesso et al.| 2006 
Loveday et al. |2012| l and redshift (e.g., |De Propris et al. | 1999j [Lin 



etal. 



2006 , Muzzin et al. 2008 , Mancone et al. 2010 Capozzi et al. 
2012 1. Since the IRAC photometry probes the peak of the stellar 
light, the IR LF serves as a good proxy for the stellar mass func- 
tion. In this section, we present the first measurements of the IR 
galaxy luminosity functions of SZ-selected clusters, and compare 
our results to those found from IR-selected cluster samples. 



The BCGs are extended objects for which the aperture cor- 
rected point source magnitudes are not good approximations to the 
total flux. We therefore use the SEXTRACTOR MAG_AUTO mag- 
nitudes, again scaling up the photometric uncertainties using the 
method of Barmby et al. (20081. Note that we do not attempt to 
deblend the BCGs beyond the level provided by SEXTRACTOR. 
While it is possible to use more sophisticated techniques to im- 
prove deblending using higher resolution imaging at other wave- 
lengths (e.g., TFIT; Laidler et al. 2007 ), we only possess relatively 
shallow ground-based optical data for the objects in our sample, 
which would require a large ^-correction to IRAC wavelengths. 



3.1 Method 

We divide the cluster sample by redshift into two equal sized sub- 
samples (0.2 < z < 0.5 and 0.5 < z < 1.1) and measure their com- 
posite luminosity functions within Rsoo- Although we have a large 
sample of spectroscopic redshifts for each cluster in addition to 
optical photometry, in most cases the optical data do not provide 
coverage out to R^qo- Therefore we use statistical background sub- 
traction to ensure consistency in our analysis across all the objects 
in the sample. This lack of optical data for some regions of most of 
the cluster fields, coupled to the low resolution of the IRAC images, 
makes star-galaxy separation difficult. Following ]Lin et aT] ( |2012| >, 
we remove the brightest stars by cross-matching the IRAC cata- 
logues with the 2MASS point source catalogue (using a 2" match- 
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Figure 3. Comparison of nij 6 found from the SZ-selected cluster sample 
used in this work with results for an IR-selected cluster sample I Mancone 
|et al.|20i0) . In both cases, the faint end slope is fixed to a = —0.8. There is 
reasonable agreement, which suggests that the galaxy populations in each 
sample are not likely to be significantly different, although 6 for the (z) = 
0.68 subsample of SZ-selected clusters is approximately 0.3 mag brighter. 
The dashed line shows the expected evolution of a solar metallicity Bruzual 
|& Char iot 1 2003 1 X — 0. 1 Gyr single burst stellar population model formed 
at Zf = 3, normalised to the |Mancone et a l. 1 2010 1 measurement at z = 
0.37. The horizontal error bars mark the redshift ranges covered by the ACT 
subsamples. 



ing radius); this should reduce the stellar contamination to around 
the 7 per cent level ( |Lin et al.|2012| >. We also remove all objects 
brighter than the BCG in each field, which further reduces stel- 
lar contamination. Remaining stars are removed statistically during 
subtraction of the background sample. Since our IRAC cluster im- 
ages do not extend to a large radius beyond /Jgoo, we use the IRAC 
source catalogue of the Extended Groth Strip (EGS; Barmby et al. 
2008 ) as the background field sample. We treat the EGS cataogue 
in the same manner as our IRAC catalogues of the cluster fields 
throughout (including, for example, removing most stellar contam- 
ination by cross matching with 2MASS). We mask out areas around 
very bright stars in some of our cluster fields, and take the reduction 
in area into account when subtracting the background contribution. 

We compute the composite luminosity functions for each sub- 
sample using the method of jColless 



corrected magnitudes (see Section 2.2 



19891. 



We use the aperture 
in the observed frame, and 
exclude the BCGs since they are not drawn from the same popula- 
tion as ordinary cluster galaxies (e.g., | Tremaine & Richstone|1977| 
|Loh & Strauss||2006| [Lin et al.||2010l We apply a modest com- 
pleteness correction (approximately 10 per cent for 21.5 < 6 < 
22.0), estimated from inserting synthetic point sources into our 
images (see Section [2~2) l. For the background sample from EGS, 
which reaches similar depth to our observations, we use Table 5 
of |Barmby et al. 1^2008} to model the completeness as a function 
of magnitude. For each cluster in each subsample, we ^-correct 
its magnitudes (assuming a t = 0.1 Gyr single-burst Bru zual &| 

3 with solar 



Chariot 



1 2003 BC03 hereafter) model, formed at Zf ; 



metallicity) to the median redshift of the subsample, and take into 
account the distance modulus between each cluster and that me- 
dian redshift. We divide the data into bins of width 0.5 mag, and 
then perform the same operations on the background galaxy sam- 
ple, before subtracting the area-scaled background contribution in 
each bin. 

We fit Schechter ( 1976) functions to the composite LFs using 



X minimisation. We fit for both the characteristic magnitude (m*) 
and the faint-end slope (a), although we also perform fits with a 
fixed to other values common in the literature (e.g., a = —0.8), in 
order to simplify the comparison with other works (since m* and a 
are degenerate). The normalisation of the Schechter function fit is 
fixed such that its integral is equal to the number of galaxies in the 
composite LF. 



3.2 Results 

Fig. [2] shows the LFs in each IRAC channel for both the low 
and high redshift subsamples. The values of m* that we derive in 
both channels for the (z) = 0.42 subsample agree at better than 
la with measurements of m* in IR-selected cluster samples by 
|Muzzin et al.| ( |2008| > and |Mancone et al.| (2010| l at similar red- 
shift (note that the faint end slope, a, was fixed to -0.8 in these 
works). We find a slightly brighter m* for the (z) = 0.68 subsample 
(nij 6 = 19.26 ± 0.07 for a = —0.8) in comparison to both Muzzin 
|et al.|p008) and |Mancone et al."| (2010), although all values are con 
sistent at better than the 2a level. If this difference is real, it may 
reflect a different timescale for the build up of the bright end of 
the luminosity function between the samples, with more massive 
galaxies being assembled at earlier times in more massive clusters 
(since we expect the average mass of the SZ-selected sample to be 
significantly larger than that of the IR-selected samples, which are 
drawn from a smaller survey area). Fig. [3] illustrates the evolution 
of ?ri 



"3.6 



in the 



Mancone et al. 



1 2010) sample in comparison to our 



results. We also see reasonable agreement between m* measured 
for our (z) = 0.68 subsample and the results of De Propris et al. 
(|2007)) at (z) = 0.75, when we set a 



-0.25, as was found in that 



work. 

The faint-end slopes that we find when fitting for both m* and 
a are consistent with that measured by |Lin et al. 1^2004) in the K- 
band from a heterogeneous sample of X-ray selected clusters at 
z = 0.1, who found a = -0.84 ±0.02, although we note that the 
uncertainties on our measurements are much larger. We conclude 
that the infrared LFs of SZ-selected clusters do not differ signifi- 
cantly from clusters selected using other methods. 



4 BRIGHTEST CLUSTER GALAXIES 
4.1 The Hubble diagram 

BCGs are the brightest galaxies in the Universe in terms of their 
stellar emission, and early work recognised their potential as stan- 
dard candles, using their Hubble diagram to estimate the deceler- 
ation parameter (e.g., Sandage 1972). More recently, the Hubble 
diagram of BCGs has been used to study their mass growth (e.g., 
|Aragon-Salamanca et ak][T998l |Burke et al.l|2000"l [Broiigh et al.| 
2002). Some semi-analytic models predict that BCGs should have 
acquired about 80 per cent of their stellar mass since z ~ 1 through 
accretion and merging (De Lucia & Blaizot 2007), although recent 
observations of BCGs in both high redshift optical and X-ray se- 
lected clusters show that most of the stellar mass in these objects 
was already assembled by z ~ 1 (e.g., Whiley et al. 2008, Brough 



Lidman et al. 



[eTaI1[2008l |Collins et al.||2009| [Stott et al.||2010| 
[20T2l >. 

Fig. [4] presents the 3.6 [im Hubble diagram of BCGs in the 
ACT SZ-selected sample, in comparison to A"-band observations 
of X-ray selected clusters (/^-corrected to 3.6 u.m, assuming a x = 
0.1 Gyr single burst, solar metallicity, BC03 model formed at zf = 
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Table 2. Properties of the BCGs; ma and m4.5 are the total apparent magnitudes (observed frame) on the AB system in the IRAC 3.6 and 4.5 urn channels 
respectively; M t is the stellar mass estimated from 1713.6 alone, assuming a T — 0. 1 Gyr burst Bruzual & Chariot (20031 model with Salpeter ( 19551 IMF and 
solar metallicity formed at zf = 3. The error bars on the stellar mass estimates correspond to assumptions of zf = 2 and Zf = 5 (systematic errors due to the 
choice of stellar population model and/or IMF are neglected). 



Cluster 


RA (J2000) 


Dec. (J2000) 


'"3.6 


'"4.5 


M„ (10 u Mq) 


ACT-CLJ0102-4915 


01 h 02 m 57?772 


-49°16'19"14 


17.204 ±0.009 


17.628 ±0.010 


16.2 


H2.8 
-3.4 


ACT-CL J0215-5212 


02 h 15 ra 12?229 


-52°12'25!'09 


17.078 ±0.006 


17.349 ±0.006 


10.9 
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Figure 4. Hubble diagram of BCGs in ACT clusters (black diamonds) in 
the 3.6 urn IRAC channel. The X^-band observations (^-corrected to the 
IRAC 3.6 urn channel assuming a T - 0.1 Gyr single burst, solar metal- 
licity, [Bruzual & Chariot 12003) model formed at z/ = 3) of two samples 
of X-ray selected clusters jEbeling et al.| 1 996] |2000| |20UT) by |Stott et al~] 
1 2008 1 are shown for comparison. The dashed line shows the expected evo- 
lution for the same Bruzual & Chariot ( 2003 1 model used to perform the 
^-correction, while the solid line shows the no-evolution line; both of these 
are normalised to match the data at z < 0. 1. 



3), taken fromlStott et al.|j2008] who measure total BCG magni- 
tudes, as in this work). One of the X-ray samples to which we com- 
pare is the Massive Cluster Survey (MACS; Ebeling et al. 



2010), which is extracted from the ROSAT All Sky Survey 



et al. 



2007 



Voges 



1999), and probes a similar mass and redshift range to the 



ACT sample. Although we have normalised the evolution tracks 
shown in this figure to match the low redshift BCG sample, it 
should be noted that the massive ACT clusters will not evolve to 
have similar masses to the low redshift clusters. 

To check if the ACT and MACS BCGs trace similar popula- 
tions, we convert our estimates of 3.6 [im BCG total magnitude to 
stellar mass, assuming a solar metallicity z = 0.1 Gyr single burst 
BC03 model, formed at z/ = 3, and do the same for the MACS 
BCGs. We compare the resulting BCG stellar mass distributions 



using the two-sample Kolmogorov-Smirnov (KS) test. When com- 
paring to the whole MACS sample, the KS test returns D = 0.44, 
with null hypothesis (that the samples are drawn from the same dis- 
tribution) probability p = 0.06. Restricting both samples to include 
only the 0.4 < z < 1.0 BCGs, the null hypothesis probability in- 
creases to p = 0.4. We conclude that the ACT and MACS BCGs 
are a similar population. 



4.2 BCG stellar mass scaling relations 

We now examine the scaling of BCG stellar mass with the prop- 
erties of the ACT clusters. We convert our estimates of 3.6 [im 
BCG total magnitudes to stellar mass, assuming a solar metallic- 
ity t = 0.1 Gyr single burst BC03 model, formed at z/ = 3. The 
uncertainty in the stellar mass estimates is dominated by the choice 
of model; since formation redshifts in the range 2 < z/ < 5 are rea- 
sonable for BCGs (e.g., |Whiley et al.|2008|[Stott et al.|2008|[20T0] ), 
we adopt the stellar mass estimates inferred from models at each 
end of this redshift range as fiducial error bars. We do not take into 
account possible systematic uncertainties, which are considerable. 
The largest systematic that can affect stellar mass estimates is the 
IMF; if we adopted a Chabrier (2003) rather than Salpeter (1955 ) 
IMF, our stellar mass estimates would be 0.24 dex lower. Simi- 
larly, adopting stellar population models with a larger contribution 
to the infrared flux by thermally pulsating AGB stars (e.g., |Maras-| 
|ton"l|2005[ |Conroy, Gunn & Whitel|2009^ would result in smaller 
stellar masses. Table [2] lists the BCG magnitudes and stellar mass 
estimates derived using our adopted BC03 model. 

Fig. [5] presents the relation between BCG stellar mass and 
cluster dynamical mass. We find mild evidence for a correla- 
tion, with Spearman rank coefficient p = 0.56, and null hypoth- 
esis (i.e. no correlation) probability p = 0.04. Using the Markov 
Chain Monte Carlo (MCMC) based method of |Kelly| ( (2007) , we 
find the slope of the relation is shallow, and poorly constrained 
(A/goo « M^ J±0A ). The intrinsic scatter in the relation at fixed BCG 
is ai ogM500 = 0.14 ±0.06. 

Many other studies have previously found a correlation be- 
tween BCG luminosity and cluster mass using X-ray and optically 
selected cluster samples (e.g., |Lin & Mohr 2004; Popessc Tet al.| 
|2007||Brough et al.|2008[|Yang et al.|2008||Mittal et al.|2009) . Like 
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Figure 5. Scaling relation between BCG stellar mass and cluster dynamical 
mass. The dashed line is a regression fit to the data using the Kelly|j2007j 
method. There is no evidence that disturbed clusters (black diamonds; de- 
fined according to the criteria described in Sifon et al. 20121 follow a dif- 
ferent trend from non-disturbed clusters (white circles). 



|Haarsma et al.|(2010) , the correlation we find is only marginally 
significant; presumably this due to the small sample size. For exam- 
ple, [Whlley£r^{2008| found a correlation between BCG A"-band 
magnitude and cluster velocity dispersion (i.e. a comparable proxy 
to the dynamical mass estimates used in this work), significant at 
the 99.9 per cent level, for a sample of 81 optically selected clusters 
in the 0.02 < z < 0.96 redshift range. We note that the wide redshift 
range covered by our sample may also be a contributing factor, if 
the relation evolves with redshift. 

Fig. [6] shows the relations between BCG stellar mass and SZ- 
signal, I^oo and vrj. These are the first measurements of these re- 
lations. In these cases we find stronger evidence for correlation 
with the BCG stellar mass. For the scaling with F500, we find 
Spearman rank coefficient p = 0.66, with null hypothesis prob- 
ability p = 0.01. Applying the Kelly (2007) regression method, 
the slope of the relation is E(z) l/i D\ y 5 oo °= M^ 2±0 6 , with in- 
trinsic scatter Oi og r 500 = 0.25 ± 0.08. For the BCG M* - y rela- 
tion, we find Spearman p = 0.74, with p = 0.002. The slope of 
the relation is E(z)~ 2 yo <* M^ 8±0,4 , and the intrinsic scatter is 
Ologyo = 0.18 ±0.05. We note that "El Gordo" (J0102-4915), the 
cluster with the largest yo an d I500 value, is a clear outlier in these 
plots. However, removing this object has no significant effect on 
the fit results. 

The extensive spectroscopic data obtained on each cluster al- 
lows their dynamical states to be classified using three different 
methods, as described in Sifon et al. (2012). A cluster is flagged 
as dynamically disturbed if it satisfies any two of (i) BCG peculiar 
velocity different from zero at more than the 2a level; (ii) BCG pro- 
jected offset from the SZ peak cluster position more than 0.2^200 ; 
and (iii) greater than 5 per cent significance level in the Dressier 
|& Shectmanj ( jl988[ l test for substructure. In both Figs. [5] and [6 
we show disturbed and non-disturbed clusters with different sym- 
bols, and see no evidence for different scaling of M500, I500, ° r >'() 
with BCG stellar mass for clusters with different dynamical states, 
though of course the sample size is small. 
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Figure 6. Scaling relations between SZ-signal and BCG stellar mass: (a) 
scaling with Y^qq; (b) scaling withyo- Self-similar evolution of the SZ-signal 
with redshift is assumed. In each panel the dashed line is a fit to the data 
using the Kelly 1 2007 1 regression method. We do not see different trends 
for disturbed (black diamonds; defined according to the criteria described 
in |Sifon et al.|2012> and non-disturbed clusters (white circles). 



Table 3. Cluster total 3.6 urn magnitudes (observed frame) and total stellar 
mass (Mj 1 q) estimated within i?5oo- The stellar mass conversion is done as 
described in Section [4~2] for the BCGs, and the same caveats noted in the 
caption to Table |2] apply. Note that any contribution to Mj' q from the ICL 
is not accounted for. 
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5 CLUSTER STELLAR MASS SCALING RELATIONS 

In this Section, we examine the scaling relations between total stel- 
lar mass within R$qq (M|qq), cluster dynamical mass, and SZ signal. 



5.1 Method 

We estimate the contribution to the total light within R500 from 
cluster galaxies other than the BCG via the procedure used to mea- 
sure the composite LFs (Section |3j. We sum the flux from the 
background-subtracted number counts for each cluster down to a 
limit of m* + 2 (where m* is the value obtained for the fit to the ap- 
propriate composite LF). We add a contribution for galaxies fainter 
than our detection limit by extrapolating the LF to m* + 5, assum- 
ing that a = —0.8 (this correction adds less than 10 per cent to the 
final stellar masses). We then convert the total flux (adding in the 
contribution from the BCG, see Section[4]l into an observed frame 
3.6 (J.m magnitude, taking into account the difference in distance 
modulus and ^-correction that were applied in estimating the com- 
posite LF. The stellar mass is then estimated assuming the same 
zj — 3 BC03 model applied in Sections [3] and [4] Finally, to obtain 
spherically deprojected values, we multiply the total stellar mass 
estimates by 0.73. This deprojection factor is obtained assuming an 

■ 2.8 (the average ob- 



NFW profile ( |Navarro et al][T997) wi th c ■■ 
tained from applying the Duffy et aTj ( [2008} c — M200 relation to 
the dynamical mass estimates), and integrating to 3.K500 along the 
line of sight. Note that our measurement does not include any con- 
tribution from the intracluster light (ICL). Table [3] lists the M|{g 
value for each cluster. 



5.2 Results 

Fig.|7]shows the scaling of M|qq with dynamical mass. While there 
is some evidence of a correlation in this plot, we find that J0237- 
4939 is a clear outlier, as it has the smallest Mf™ in the sample but 
is the second ranked cluster in the sample in terms of its dynamical 
mass. If this cluster is excluded, we obtain Spearman rank coef- 
ficient p = 0.76 with null hypothesis probability p = 0.002. Ap- 
plying the Kelly (20071 regression method with J0237 excluded, 
we find the relation M500 x M^ 09 * 0,4 , with intrinsic scatter 
°logAf 5 oo = 0.10 ±0.06. If J0237 is included in the sample, the cor- 
relation is not significant (p = 0.47, p = 0.09) and a much shallower 
slope for the relation is inferred (M500 1 



Mstar0.3±0.4\ 



'500 

We present the relations between Mf^ and SZ-signal (yo, 
T500) in Fig. [8] In both cases we see a correlation, although the 
Spearman rank test indicates that the correlation with F500 is more 
significant (p = 0.63, p = 0.02) compared to the correlation with 
>0 (P = 0.46, p = 0.10). In both cases, as with the scaling with the 
BCG stellar mass, the slopes of the relations are poorly constrained 
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Figure 7. Scaling relation between total stellar mass (M^) and dynamical 
cluster mass (Msqq). The dashed line indicates the fit obtained with the | Kelly] 
( 2007 1 regression method, after excluding the outlier J0237-4939, which 
has a very low stellar mass compared to its dynamical mass. 
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6 DISCUSSION 

We now discuss the variation of the baryon fraction in stars (/f Q = 
M|qq/M5oo) with cluster mass for the ACT SZ-selected sample. 
This relationship may give some insight into the role of feedback 
in regulating star formation efficiency in clusters (e.g., Bo de et al.| 
|2009l >. 

Fig. [9] shows our results in comparison to a number of works 
in the literature. For the ACT sample, we find /™ values that span 



(M. 



Figure 8. Scaling relations between SZ-signal and cluster stellar mass: (a) 
scaling with F500; (b) scaling with yo- Self-similar evolution is assumed in 
both cases. The dashed line indicates the fit obtained using the |Kelly|{2007) 
regression method. 
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Figure 9. Fraction of the total cluster mass in stars (/Jm) or total baryon fraction (J^qq, '- e -< ^500 P' us § as mass > where available) within R500, as a function 
of cluster mass. The solid line marks the cosmic baryon fraction iKomatsu et al. 2011 the shaded area indicates the uncertainty). The legend indicates 



Lin et al. 



2012 



are scaled up from a |Kroupa 



2001 



IMF 



measurements of /BS 

by adding 0. 1 3 dex to the total stellar masses. Only the work of Gonzalez et al. 1 2007 1 includes the contribution from the ICL. The five black diamonds indicate 
/ s b 00 estimates for J0102-4915 (gas mass measurement taken from |Menanteau et al.|2012a , J0509-5341, J0528-5259, J0546-5345, and J0559-5249 (gas mass 
measurements taken from Andersson et al. 2011} . 



the range 0.006-0.034, with median 0.022. The cluster at the low 
end of this range is J0237-4939, which has an unusually low stel- 
lar mass given its total mass. It is flagged as a disturbed cluster 
according to the criteria used by Sifon et al. (2012), although the 
other disturbed clusters in the sample do not exhibit such low 
values. In Fig.|9]we see the decreasing trend of with increasing 
cluster mass seen in other works, though the mass range covered by 
the ACT sample is relatively narrow. The values we find are 
typically slightly larger at a given mass in comparison to the results 
of other studies shown in the plot. While this could be a real dif- 
ference, there are several possible sources of systematic error, or 
differences in the analyses or assumptions applied in each study, 
which could account for this. 

The single largest source of systematic error is the choice of 
IMF. Here, we have adopted the |Salpeter|fl955| > IMF (as did the 
work of |Giodini et al.||2009| l, and in Fig. 9] we have rescaled the 
results oflLirTeT^LTi 2012} from thejKroupa (2001) IMF assumed 



in that work to Salpeter by adding 0.13 dex to their M™ values. If 
we had instead adopted a Chabrier ( 2003 ) IMF instead of Salpeter 
our stellar mass measurements would be 0.24 dex lower. Gonzalez 
|et al.| ( |2007[ > took a different approach, converting light to stellar 
mass using an empirically determined sub-Salpeter M/L ratio taken 
from observations of local ellipticals and SOs by Cappellari et al. 
(2006 ). We note that the possible variation of the IMF according 
to galaxy type is currently under debate; some recent studies jvan| 
Dokk um & Conroy|2010)|Cappellari et al.|2012") suggest that early 
type galaxies have a more bottom-heavy IMF relative to late type 
galaxies. Accounting for such an effect is beyond the scope of this 
paper. 

In a similar way, the assumption made in converting light to 
mass is also a potential source of systematic error. In the absence of 



multi-band photometry covering the entire sample out to ^500, we 
have simply assumed a constant mass-to-light ratio for each clus- 
ter, assuming all member galaxies to be represented by a single- 
burst BC03 model with z/ = 3. This is motivated by the fact that 
early type galaxies form the dominant population in massive clus- 
ters, even beyond z > 1, with the scatter about the red-sequence 
in such clusters being consistent with the bulk of star formation 
having taken place at z ~ 3 (e.g., |Mei et a l. 2009, Hilton et al. 



|2009| |Strazzullo et al. 2010). This assumption is also made in |Lin 
|et al.|p012} , and there is reasonable agreement between the two 
studies in the mass range of overlap. However, some other studies 
take into account the variation in mass-to-light ratio with galaxy 
ty pe. |Giodini et aL] ([2 009 ) apply the A^-to-stellar mass relations 



1 2007 1, which are derived from SED fitting to 



of |Arnouts et al, 
multiband photometry, according to galaxy colour. Leautha ud et al.| 
(2012 ) used SED fitting to COSMOS photometry to obtain stellar 
mass estimates for each individual galaxy, and estimated the stel- 
lar fraction using a statistical halo occupation distribution model 
approach, finding much lower values (by a factor of 2-5) than for 
other studies, including ours (although much of this difference can 
be attributed to the adoption of the |Chabrier|^2003| ) IMF). We note 
that such a galaxy-type dependent estimation of stellar mass will 
have a much larger impact at the group scale, the focus of the work 
by |Giodini efaL] ( |2009) and |Leauthaud et al.|{20T2| >, in comparison 
to the massive clusters we consider here. 



As noted earlier, we do not account for the presence of ICL 
when computing M^qq. Of the studies shown in Fig. |5] only the 
work of Gonzalez et al. ( 2007 ) measured and included the ICL 
component, which was found to be approximately 30 per cent of 
the total light across a sample of z ~ 0. 1 clusters. Recent studies 
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suggest that the ICL makes up a smaller fraction of the total cluster 
light at high redshift 4 per cent at z = 0.8; |Burke et al.|2012| >. 

Despite the limitations of the analysis in this work, the stel- 
lar fractions we find are nevertheless in reasonable agreement with 
the results of numerical simulations by Batta glia et al.| ( |2012| >. 
These cosmological simulations include sub-grid models for radia- 
tive cooling, star formation, and AGN feedback (Battagli a et al.| 
|2010j >. As can be seen in Fig. [9] our observations (in common with 
the other observational results plotted) show a larger scatter than the 
simulations, which suggests that the sub-grid models in the simula- 
tions need to be improved to capture the larger observed variations 
in/|$(M 500 ). 

The stellar component of clusters makes only a small con- 
tribution to the total baryonic mass. Measurements of gas mass 
within i?5oo (Affoo) are available from the literature for five of 
the objects in the ACT sample (Andersson et al. 2011, Menan- 
|teau et aT]|2012a) . In Fig. [9] we plot the total baryon fractions 
C/500 = [^500 +^foo]/^50o) for these clusters, in comparison to 
the cosmic mean as measured by WMAP (Ko matsu et al.||201 1| >. 
We see no evidence for a shortfall of baryons in these objects given 
current uncertainties, although this is clearly dependent upon the 
validity of the assumptions made in estimating /j™ . 



7 CONCLUSIONS 

We have performed the first study of the stellar mass component of 
galaxy clusters selected via the Sunyaev-Zel'dovich effect, using 
Spitzer Space Telescope infrared observations of a sample detected 
by the Atacama Cosmology Telescope. We found: 

(i) The 3.6 [im and 4.5 vim luminosity functions are similar to 
those measured for IR-selected cluster samples. We measure both 
a characteristic magnitude (in*) and faint-end slope (a) for the 
low redshift (0.2 < z < 0.5) ACT sample similar to those found 
by Mancon e et al.| J2010) , while for the higher redshift sample 
(0.5 < z < 1.1), the ACT clusters have slightly brighter m* (at the 
2a level). 

(ii) The relation between BCG stellar mass and cluster dynam- 
ical mass for the ACT sample is shallow (M 500 °= M° 7±0 ' 4 ). The 
relation is not very well constrained, most likely due to the small 
sample size. There is strong evidence for the correlation of BCG 
stellar mass with SZ-observables, although again the constraints 
we obtained on the slopes of these relations are poor. For the scal- 
ing with Y500, we found E(z) ~ 2 / 3 T)\ K500 °= M' 2±0 ' 6 , with intrinsic 
scatter O iogY500 = 0.25 ± 0.08. 

(iii) Excluding the cluster J0237-4939, which has anomalously 
low total stellar mass (M|qq) compared to its dynamical mass, we 
found M500 °= Mf^°- 9±0A , with intrinsic scatter Gi og M 5m = 0.10± 
0.06. We also made the first measurement of the scaling of Mi^ 
with SZ-signal, finding £(z)~ 2 / 3 D\ Fg o Mfgg 10±0 - 6 . 

(iv) The stellar fractions that we measured cover the range 
0.006-0.034, with median 0.022. These are larger than found in 
some other studies of clusters of similar mass, but in reasonable 
agreement with the results of a similar analysis of X-ray selected 
clusters using WISE data ( |Lin et aL| [2012 1. For the five clusters 
with additional gas mass measurements available in the literature, 
we see no evidence for a shortfall of baryons in clusters relative to 
the cosmic mean value. 

In the future, we intend to extend this study to include ACT 
clusters on the celestial equator (Hasselfield et al. 2012; Me nanteau] 



|et al.|2012b[ >, where the overlap with deep SDSS Stripe 82 optical 
photometry ( |Ann is et al. 2011) will allow us to perform a more 
detailed analysis, without the need for some of the assumptions 
used in this work. 
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